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Abstract 

We review recent lattice computations relevant for D and B decays. De- 
cay constants fo,D a , Ib,b 3 , D — > K(K*), B — > tt, p semi-leptonic form 
factors together with the slope of the Isgur-Wise function calculations are 
presented. Some recent results of B — > K*j form factors will be given. 
1/M corrections to the asymptotic scaling laws are discussed. 

1 Introduction 

B meson physics is now an active field (where several thousand of physi- 
cists are engaged) in which lattice QCD takes advantage from recent im- 
provements increasing the statistics and improving the control of the sys- 
tematics through theoretical progress (improved actions- • ■). In addition 
to phenomenological predictions, lattice QCD can test the scaling laws 
predicted by HQET (Heavy Quark Effective Theory). Up to now, precise 
predictions can be made only in the quenched approximation for which 
the systematic errors cannot be evaluated; however, results obtained in 
the unquenched theoryjlj suggest that the quenching effect is small when 
one deals with heavy quarks. In general, comparison with experiment 
shows an agreement with lattice data inside error bars which are still siz- 
able although decreasing. This is a review of recent calculations done in 
the quenched approximation using the Wilson and clover (continuum limit 
improved) actions on heavy-light meson decays. 



2 Strategy to study B meson on the lattice 

Since the inverse lattice spacing a -1 ranges from 2 to 4 GeV, one can- 
not study B meson directly on the lattice. Indirect informations can be 
obtained through the following: 

- One uses a set of relatively heavy mesons with masses up to 0.7a -1 , 
i.e. heavier than the D but lighter than the B meson (" fictitious D 
mesons"). We call this mass region the "moving quark" region. 
- On the other hand, a method proposed by Eichtenj| allows one to put 
a quark with infinite mass on the lattice and the latter is considered in 
this approach as a static source of color. We call this mass region the 
"static quark" one (where heavy flavours are studied at lowest order in jj 
expansion, M being the heavy meson mass). 
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A physical quantity computed in these two mass regions is interpolated 
to the B meson with the help of the scaling laws of the HQET. The value 
in the static limit reduces the uncertainties due to the extrapolation. This 
method has shown to be very effective (Fig.l) in the estimation of the 
decay constants. In semi-leptonic decays, the calculation in the static 
limit is not yet available, but one can study the scaling behaviour and try 
an extrapolation to the B meson. The predictions concerning the B meson 
semi-leptonic decays remain at a semi-quantitative level but it shows that 
the extrapolation may be done and improvements in the near future are 
expected. The next section is an illustration of the method explained 
above. 



3 Leptonic decays D(B) — > tv 



A hadron mass can be obtained from the study of an appropriate Eu- 
clidean correlation function as the coefficient of its exponential time de- 
pendence: 

G{t) = /d 3 x(u(x,i) 7o75 c(x,i)c(0, 0) 7 o7 5 u(O,0)) 

l e -m D t ^ 



fr> m D „- mD / 



The determination of the expectation value in eqjl] is a non perturbative 
problem which can be solved numerically. The second approach (static) 
is based on the expansion of the heavy quark (7?) propagator in inverse 
powers of the quark mass as proposed by EichtenQ; the H is static and 
does not live effectively on the lattice but the quantity f H -J M H can be 
measured and is predicted to be independent of the heavy mass. The 
confrontation between the two methods is presented in Fig.l. The HQET 
tells us that when Mjj — > oo, the vector (V) and pseudoscalar (P) decay 
constants scale with the mass of the heavy quark, Mh, 0-H(M = Mp — 
My = Mjj, A) = 11 - f JV» as: 

^=f P = ^=a 8 {M)-^. (2) 
fv VM 
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Ref. 


(3 


MMeV) 


/c.(MeV) 


ELC(W) 




6.4 


210 ±15 


227 ± 15 


APE (C) 




6.0 


218 ±9 


240 ±9 


BLS(W)[ 


'I 


6.3 


208(9) ±35 ±12 


230(7) ±30 ± 18 


UKQCD (C)g 


6.2 


185+^ +f 


212™ b 


WA75|9| 






232 ± 45 ± 20 ± 48 


CLE02|n| 






344 ± 37 ± 52 ± 42 


ARGUS |jj 






267 ± 28 



Table 1: W (C) refers to Wilson (clover) action. 
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Figure. 1. Linear and quadratic fit in 1/M are shown, the vertical line shows the physical 
B meson. 

In Figure. 1, we notice the consistency between the moving quark results 
and the static ones. It appears that there are large corrections to the 
asymptotic scaling behaviour (2). The lattice results compared to the 
experimental ones are reported in Table 1(2) concerning the D(B) meson. 

From Table 1, one can see that the different lattices agree more or less. 
Up to 5% - 10% we find: f D ~ 210MeV and f Ds ~ 230MeV. f Ds which 
has been predicted by lattice since several years, will provide an important 
check since the large experimental errors may be substantially reduced in 
the future. 

In Table 2, the general tendency is f B ~ 200MeV (up to 20%), in 
agreement with QCD Sum Rules calculations. 

The B— parameter: the predictions for the B — B mixing depend on the 
B— parameter of the heavy light AB — 2 four quark operator. ELC 

givesg): B d0 = 1.05 ± 0.08, B g0 = 1.16 ± 0.07 and ~ = 

D d B d 

1.02±0.02; the prediction for the physically relevant combination in B — B 
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Ref. 


P 




Jb s 


ELC(W)|4] 


6.4 


205 ± 40 


1.08 ± 0.06 


APE (S-C)|12[ 


6.2 


290 ± 15 ±45 


1.11(3) 


APE (S)[14l 


6.0 


350 ± 40 ± 30 


1.14(4) 


APE (S-C)|4 




6.0 


328 ± 36 


1.19(5) 


UKQCD (C)[S 


1 


6.2 


iou_ 6 _ 19 


1.22 ±0.04 


UKQCD (S-C) 


i 


6.2 


z00 -15 -14 


1.141* 


BLS (S)|7I_ 


6.3 


235(20) ± 21 


1.11 ± 0.05 


Allton(S-W)[l| 


6.0 


310 ±25 ±50 


1.09 ±0.04 


HEMCGCfl] 


5.6 


200 ± 48 





Table 2: S refers to Static limit. The HEMCGC unquenched result agree with 
those obtained in the quenched approximation. 



mixing and CP violation is: f Bd y / B Bd = 220 ± 40MeV. 

4 Semi-leptonic decays D(B) —> K, K*(ir, p)iv 

In this study, we first calculate D — > K, K* then extrapolate to B — > 
7r, p. The amplitudes are expressed in terms of four form factors f + , V 
and A\ i and we need to know them at momentum transfer q 2 = (where 
we have the maximum of phase space). The present lattices run in the 
close vicinity of q 2 = for the D meson, far from so for the B one. In the 
latter case we proceed in two steps: 

- The first is to extrapolate to the B meson mass at fixed momentum 
near the no recoil point (g^ax) anc ^ this is doable with the help of HQET: 
when M — > oo at fixed q and \\q\\ <§; M, the form factors scale as|]i"6||: 
/+, V, A 2 ~ M 1 / 2 , and A 1 - M" 1 / 2 ). 

- The 2 nd and difficult step is to extrapolate to q 2 = 0; people often use 
the Nearest Pole Dominance approximation VMD (F(q 2 ) — j^pJW 1 ^ ^* ^ s 
the exchanged meson mass in the t-channel) which has no firm theoretical 
grounding. 

D meson study: lattice results compared to model predictions and exper- 
imental data are reported in Table 3; the two first lines |ll| , |L9| have been 
recently obtained using the clover action. There is agreement between 
lattices, models and experiments for / + and A\. For V the central value 
is below within errors. Concerning Ai, there used to be a problem but 
now the central value agrees with the experimental average although with 
large errors. 

B meson study : ELC(T| and APE@ have studied the HQET scaling 
laws near q 2 lax and it seems that there are large corrections for V and 
Ai t 2- The B meson form factors predictions at q 2 = have large errors 
and rely on the nearest pole dominance assumption p7[ , |l8|] . This type of 
extrapolation has been done for the first time and now, one knows how to 
handle -h corrections. 

The Isgur-Wise function £(x): when the final meson is heavy, semi- 
leptonic form factors are expressed in terms of one universal function, 
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Ref. 


f+(0) 


V(0) 


Ax(0) 


EXP pi 


0.77(4) 


1.16 ± .16 


0.61(5) 


APE[ 

ukqc: 

ELC 
APE 
BLS[ 


18| 

13 
11 

20f 


0.72(9) 

U.O 1 _g 

0.60 ±0.15 ±0.07 

0.63 ±0.08 
0.90 ±0.08 ±0.21 


1 AH I A OA 

1.00 ± 0.20 

0.86 ±0.24 
0.86 ±0.10 
1.43 ±0.45 ±0.49 


0.64 ±0.11 

u -' u -io 
0.64 ±0.16 

0.53 ±0.03 

0.83 ±0.14 ±0.28 


SR|23] 


60 +u ib 
u - ou -o.io 


1.10 ±0.25 


0.50 ±0.15 


QM.l 
QM.2 


21 

22 




U. (0 

U.o 


1.26 
1.1 


A OO 

U.OO 
U.o 


Ref. 


A 2 (0) 


V(0)/A!(0) 


A 2 (0)/Ax(0) 


EXP§J 


0.45(9) 


1.90 ±0.25 


0.74 ±0.15 


APE[ 
UKQC] 
ELC 
APE 
BLS[ 


17 
15 
20 


LS] 


0.46 ±0.34 

0.40 ±0.28 ±0.04 

0.19 ±0.21 
0.59 ±0.14 ±0.24 


1.59 ±0.29 

1.3 ±0.2 
1.6 ±0.2 
1.99 ±0.22 ±0.33 


0.73 ±0.45 

0.6 ±0.3 
0.4 ±0.4 
0.7±0.16±0.17 


SR|23f 


0.60 ±0.15 


2.2 ±0.2 


1.2 ±0.2 


QM.l 
QM.2 


21 

22 




1.15 
0.8 


1.4 
1.4 


1.3 
1.0 



Table 3: D — > K,K* semi-leptonic form factors. EXP, LAT, QM and SR refer 
to experimental average, lattice, quark model and QCD sum rules calculations 
respectively. 

unknown apart from the no recoil point £(1) = 1. In practice we try to 
measure its slope around this point p 2 = — £'(1). In the D mass region, 
UKQCD found p 2 = 1.2^|| while CLE02 gives p 2 = 1.01 ± 0.15 ± 
0.09||. A preliminary analysis (UKQCD-MartineUi0) gives p 2 = 1.7 ± 
0.2 for masses m > mo,- 

5 Radiative Decay B — > K*^y 

The decay rate is expressed in terms of Ti 2 (q 2 ) where for a real photon 
(q 2 = 0), we have the exact condition 7i(0) = T 2 (0); the method used to 
extract Ti. 2 is the same as the semideptonic one: 

- Use the HQET scaling rules near q\ ax (at leading order, Ti~\/M, 
T 2 ~-^) to extrapolate to M B . With the clover action, APE(/3 = 6.0) ||§ 

and UKQCD(/3 = 6.2)|§ find for T 2 {q 2 lax ): 0.21(2) and 0.269±g 7 respec- 
tively. 

- The problem is: how to extrapolate to q 2 = 0? If we apply VMD 
on the two factors simultaneously, we find that they scale differently at 
q 2 = (Ti(0) ~ M- 1 ' 2 and T 2 (0) ~ A/~ 3 / 2 ) in contradiction with the 
fact that they must be equal at q 2 = 0. So for which one is VMD better? 

* Applying VMD on T 2 , APEp|, UKQCD|9| and Bernard et al.,@ 
find for T 2 (0): 0.084(7), 0.112l^^and 0.10 ± 0.01 ± 0.03 respectively. 

* Applying VMD on T u APE@ finds a larger value Ti(0) = 0.20(7). 

There is thus a contradiction between the two approaches; it seems 
to us that the higher value is favoured because indications from lattices 
(APE and UKQCD) show that the q 2 dependence of T 2 is much weaker 
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than would be predicted by VMD. Using the 1 st approach, UKQCD finds 
BR(B -> K* 7 ) = (1.7±0.6(stat)l^(sys) x 10~ 5 while APE, when applying 
VDM on Ti finds a value closer (preliminary) to CLEOj|l| result: BR = 
(4.5 ± 1.5 ±0.9) x 10~ 5 . 

6 Conclusion 

There have been good quantitative studies for fxj s (~ 230MeV), fr> (~ 
210MeV), f B (~ 200MeV (20%)) and for D -> K^lv (agreement up to 
10%). The study of jj corrections (which are found to be large) to HQET 
are now under control. Concerning the B — > Tt(p)£is and B — > K*j, since 
the q 2 behaviour is still largely unknown, the predictions at q 2 = are 
only qualitative. 
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